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§ : Abstract 

A moving frame formulation of non-stretching geometric curve flows in Euclidean space is used 
to derive a 1+1 dimensional hierarchy of integrable S'0(3)-invariant vector models containing the 
Heisenberg ferromagnetic spin model as well as a model given by a spin-vector version of the 
mKdV equation. These models describe a geometric realization of the NLS hierarchy of soliton 
equations whose bi-Hamiltonian structure is shown to be encoded in the Frenet equations of the 



moving frame. This derivation yields an explicit bi-Hamiltonian structure, recursion operator, and 
constants of motion for each model in the hierarchy. A generalization of these results to geometric 



CO 



surface flows is presented, where the surfaces are non-stretching in one direction while stretching in 



all transverse directions. Through the Frenet equations of a moving frame, such surface flows are 

o ; 

shown to encode a hierarchy of 2+1 dimensional integrable S*0(3)-invariant vector models, along 

with their bi-Hamiltonian structure, recursion operator, and constants of motion, describing a 
•i— i , 

geometric realization of 2+1 dimensional bi-Hamiltonian NLS and mKdV soliton equations. Based 

5-H ' 



on the well-known equivalence between the Heisenberg model and the Schrodinger map equation 
in 1+1 dimensions, a geometrical formulation of these hierarchies of 1+1 and 2+1 vector models 
is given in terms of dynamical maps into the 2-sphere. In particular, this formulation yields a new 
integrable generalization of the Schrodinger map equation in 2+1 dimensions as well as a mKdV 
analog of this map equation corresponding to the mKdV spin model in 1+1 and 2+1 dimensions. 
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I. INTRODUCTION AND SUMMARY 



Spin systems are an important class of dynamical vector models from both physical 
and mathematical points of view. In physics such models describe nonlinear dynamics of 
magnetic materials, while in mathematics they give rise to associated geometric flows of 
curves where the unit tangent vector along a curve is identified with a dynamical spin 
vector. 

A main example [14] is the Heisenberg model for the dynamics of an isotropic ferromagnet 
spin system in 1+1 dimensions. The geometric curve flow described by this S'0(3)-invariant 
model corresponds to the equations of motion of a non-stretching vortex filament in Eu- 
clidean space. Remarkably, the vortex filament equations are an integrable Hamiltonian 
system that is equivalent to the 1+1 dimensional focusing nonlinear Schrodinger equation 
(NLS) through a change of dynamical variables known as a Hasimoto transformation [13]. 

The vortex filament equations are one example in an infinite hierarchy of non-stretching 
geometric flows of space curves whose equations of motion have a well-understood integra- 
bility: e.g. a Lax pair and an associated isospectral linear eigenvalue problem; an infinite 
set of symmetries and constants of motion; and exact solutions with solitonic properties. 
This integrability structure turns out to have a simple geometric origin. In particular, all of 
these equations of motion are generated through a recursion operator that can be derived 
geometrically [10, 23] from the Serret-Frenet structure equations given by a 5*0(3) moving 
frame formulation for arbitrary non-stretching curve flows in Euclidean space, with the com- 
ponents of the frame connection matrix providing the dynamical variables that appear in the 
equations of motion. More recently, these SO (3) frame structure equations have been found 
to geometrically encode a pair of compatible Hamiltonian operators that yield a concrete 
bi-Hamiltonian structure for the equations of motion of each integrable curve flow in the 
hierarchy [17]. 

The explicit bi-Hamiltonian form of the resulting equations of motion depends on a choice 
of the SO (3) moving frame for the underlying space curve, which determines the form of the 
frame connection matrix and hence yields the dynamical variables in terms of the curve. In 
the case of the vortex filament equations, the dynamical variables consist of the curvature 
invariant, k,, and the torsion invariant, r, of the space curve, corresponding to the choice of 
a classical Frenet frame [12] given by the unit tangent vector, unit normal and bi-normal 
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vectors, along the curve. Other geometrical choices of a moving frame can be made [2], 
since there is a 5*0(3) gauge freedom relating any two orthonormal frames along an arbitrary 
curve in Euclidean space. In particular, the Hasimoto transformation arises geometrically as 
a gauge transformation from a Frenet frame to a parallel frame [7] , where the frame vectors 
in the normal space of the curve are chosen such that their derivative with respect to the 
arclength s along the curve lies in the tangent space of the curve. This choice of frame is 
unique up to rigid SO (2) rotations acting on the normal vectors by the same angle at all 
points along the curve, while leaving invariant the tangent vector. The corresponding pair 
of dynamical variables (defined by the connection matrix of a parallel frame) are naturally 
equivalent to a single complex-valued variable u = K,exp(i J rds) that is determined by the 
curve only up to constant phase rotations u — > e l ^u (where is independent of arclength s). 
This dynamical variable u thus has the geometrical meaning [5] of a U (1) ~ SO (2) covariant 
of the space curve. Importantly, the resulting Hamiltonian structure for the equations of 
motion looks simplest in terms of the covariant u, which directly incorporates the Hasimoto 
transformation, rather than using the classical invariants k and r. 

The purpose of the present paper will be to give some new applications of these ideas to 
the study of integrable vector models in 1+1 and 2+1 dimensions. 

Firstly, from the hierarchy of non-stretching geometric space curve flows that contains 
the vortex filament equations, we derive the complete hierarchy of corresponding integrable 
S'0(3)-invariant vector models in 1+1 dimensions, along with their bi-Hamiltonian integra- 
bility structure in explicit form. In addition to the Heisenberg model, this hierarchy will 
be seen to contain a model that describes a spin-vector version of the mKdV equation. 
Our results provide a new derivation of the Hamiltonian structure, recursion operator, and 
constants of motion for these models. 

Secondly, we extend the derivation to a geometrically analogous class of surface flows 
where the surface is non-stretching in one coordinate direction while stretching in all trans- 
verse directions. Such surfaces arise in a natural fashion from a spatial Hamiltonian flow 
of non-stretching space curves. This generalization will be shown to give rise to a class of 
2+1 dimensional NLS and mKdV soliton equations with an explicit bi-Hamiltonian struc- 
ture, yielding a hierarchy of integrable SO (3)-invariant vector models in 2+1 dimensions. 
In particular, this hierarchy includes 2+1 generalizations of the Heisenberg spin model and 
the mKdV spin model, which were found in earlier work by one of us [15, 18-21]. Our 
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derivation here, in contrast, yields the explicit bi-Hamiltonian structure, recursion operator, 
and constants of motion, which are new results for these models. We also write out the cor- 
responding surface flows explicitly in terms of geometric variables given by [12] the geodesic 
and normal curvatures and the relative torsion of the non-stretching coordinate lines on the 
surface. The surface flow arising from the 2+1 integrable Heisenberg model will be seen to 
describe a sheet of non-stretching vortex filaments in Euclidean space. 

Lastly, we also derive an interesting geometric formulation of these results by viewing 
the spin vector as a dynamical map into the 2-sphere in Euclidean space. This formula- 
tion is based on the well-known geometrical equivalence between the Heisenberg model and 
the Schrodinger map equation in 1+1 dimensions [25]. When applied to the 1+1 and 2+1 
dimensional hierarchies of S'0(3)-invariant vector models, our derivation yields a new inte- 
grable generalization of the Schrodinger map equation in 2+1 dimensions as well as a new 
mKdV analog of this map equation corresponding to the mKdV spin- vector model in 1+1 
and 2+1 dimensions. 

The rest of the paper is organized as follows. In section II, we review from a unified point 
of view the mathematical relationships amongst 1+1 dimensional vector models, dynamical 
maps into the 2-sphere, non-stretching curve flows in Euclidean space, Frenet and parallel 
frames, and the Hasimoto transformation. In section III, we derive the NLS hierarchy of 
soliton equations in terms of the geometrical covariant u given by the Frenet equations of 
a moving parallel frame for non-stretching space curve flows. This approach directly yields 
the explicit bi-Hamiltonian structure of these soliton equations, including a formula for the 
Hamiltonians. As examples, the parallel-frame Frenet equations are used to show, firstly, 
how the NLS equation itself corresponds geometrically to the Heisenberg spin model and 
the Schrodinger map equation; and secondly, how the mKdV spin model and the mKdV 
map equation arise geometrically from the next soliton equation in the NLS hierarchy. 

Section IV contains several main results. We work out the equations of motion for the 
space curves corresponding to the NLS hierarchy and write down the induced flows on the 
curvature and torsion invariants k, t. Next we derive the resulting geometrical hierarchies 
of vector models and dynamical map equations, along with their bi-Hamiltonian structure, 
recursion operators, and constants of motion. This new derivation involves only the parallel- 
frame Frenet equations plus the bi-Hamiltonian structure of the NLS hierarchy. The explicit 
bi-Hamiltonian form of the Schrodinger map equation and Heisenberg model, including a 
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geometric expression for the Hamiltonians, are presented as examples. 

In section V, we consider surfaces generated by a spatial Hamiltonian flow of curves with a 
parallel framing in Euclidean space. The underlying Hamiltonian structure is shown to arise 
naturally from the Frenet equations of the induced frame along the surface. This formulation 
is then used in section VI to study surface flows expressed in terms of the covariant variable u 
geometrically associated with the non-stretching space curves that foliate the surface, where 
the surface is stretching in all directions transverse to these curves. We show that the bi- 
Hamiltonian structure for 1+1 flows on u has a natural extension to 2+1 flows based on the 
observation that the Hamiltonian operators involve only the coordinate in the non-stretching 
direction on the surface. This leads to a hierarchy of 2+1 flows on u, with the starting flow 
given geometrically by translations in the coordinate in the transverse direction, which yields 
a 2+1 generalization of the NLS hierarchy. 

The final two sections of the paper contain our main new results. In section VII, we 
use the surface Frenet equations to derive the complete hierarchies of integrable 2+1 vector 
models and dynamical maps arising from the 2+1 generalization of the NLS hierarchy. The 
derivation yields the explicit bi-Hamiltonian structure of these two hierarchies, in addition 
to their respective recursion operators and constants of motion. As examples, the integrable 
generalizations of the Heisenberg model and the mKdV spin model in 2+1 dimensions are 
written down in detail, as well as the corresponding new 2+1 dimensional integrable gen- 
eralizations of the Schrodinger map equation and mKdV map equation. In section VIII, 
we work out the equations of motion for the surface flows that correspond to the previ- 
ous hierarchies. These equations are obtained by means of a different framing defined in 
a purely geometrical fashion by the non-stretching coordinate direction on the surface and 
the orthogonal direction of the surface normal in Euclidean space. We also discuss aspects 
of both the intrinsic and extrinsic geometry of the resulting surface motions. In particular, 
we obtain a recursion operator, constants of motion, and explicit evolution equations for- 
mulated in terms of geometric variables given by the geodesic curvature, normal curvature, 
and relative torsion of the non-stretching coordinate lines on the surface. 

Some concluding remarks on future extensions of this work are given in section IX. 
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II. VECTOR MODELS AND SPACE CURVE FLOWS 



We start from an arbitrary 5*0(3) vector model in 1+1 dimensions, 

S t = f(S,S x ,S xx ...), \S\ = 1 (2.1) 

where S(t,x) = (Si, £2, S3) is a dynamical unit vector in Euclidean space, / is a vector 
function _L S, and x belongs to some one-dimensional domain C. A running example will 
be the Heisenberg spin model 

S t = S A S xx = (S A S x ) x (2.2) 

with C being 1 or S 1 . 

There are two different ways to associate a curve flow to equation (2.1). One formulation 
consists of intrinsically identifying S with a map 7 into the unit sphere S 2 C K 3 . Then St 
and S x correspond to 74 and 7^; d x + S(S X -) corresponds to the covariant derivative V x on 
the sphere with respect to the tangent direction 7^; and SA corresponds to the Hodge dual 
* = J (i.e. a complex structure on the sphere). Under these identifications, each vector 
model (2.1) describes a curve flow 

lt = F( lx ,V xlx ,...) (2.3) 

for j(t,x) on S 2 . Ex. the Heisenberg model (2.2) corresponds to 

It = JVxlx (2.4) 

which is the Schrodinger map equation on S 2 . 

Alternatively, in an extrinsic formulation, S can be identified with the unit tangent vector 
T along a non-stretching space curve given by a position vector f in Euclidean space, 

S = T = r x (2.5) 

where x is the arclength along the curve f(x). Then the equation of motion of r is 

r tx = f(f x ,f xx , . . .), I f x I = 1, (2.6) 

or equivalently 

/X 
f(f x ,f xx ,...)dx, \f x \ = 1, (2.7) 



under which the arclength of the curve is preserved, i.e. f c \r x \dx = £ is a constant of the 
motion. Ex. the Heisenberg model (2.2) corresponds to 



r t = r x A r x 



(2.8) 



with \ f x \ = 1. This is the equation of motion of a non-stretching vortex filament studied by 
Hasimoto [13]. 

To proceed we first introduce a Frenet frame E along f(x). It is expressed in matrix 
column notation by 

M 



E = 



N 



(2.9) 



where 



T t x , N | T x | T x | t xx | t xx , B T A A/ - | t xx \ T' x A r x:r . 



(2.10) 



Here iV is the unit normal and B is the unit bi- normal of the space curve r(x). Note we 
have the relations 

S = T, S x = kN, SAS x = kB, (2.11) 

where 



k — T r • N — \S T 



(2.12) 



is the curvature of f(x), and 

r = N x -B= \S X \~ 2 S XX ■ (S A S x ) (2.13) 
is the torsion of f(x). The Serret-Frenet equations of this frame (2.9) are given by 

E,, = KE (2.14) 

with 



K 



\ 



Gso(3). 



k 

—K T 

-r 

From the equation of motion (2.1) for S we obtain the frame evolution equation 



V 



(2.15) 



/ 



a 2 a 3 ^ 



E, = AE, A 



-a 2 ai 
-a 3 — a\ 



Gso(3) 



(2.16) 



where 



a 1 = f x -B/\S x \ = f x -(SAS x )/\S x \ 2 , 

a 2 = f-N = f-S x /\S x \, 

a 3 = f-B = f-(SAS x )/\S x \, 

are determined by taking the t-derivative of (2.11) and substituting (2.1), followed by ap- 
plying respective projections orthogonal to T, N, B. 

This evolution of the frame E induces evolution equations for k and r through the zero- 
curvature relation K t = A x + [A, K]. Ex. the Heisenberg model (2.2) gives the vortex 
filament equations in terms of the curvature and torsion [14]: 

9 _ (^t) x 

Kt — KT X ^H X T — , 

K 

f^xxx K> XX K X / K m 2 i 1 2\ 

n = 5 2tt x + kk x = ( r + -k ) x . 

K K K 2 

(2.17) 

Next we perform a 5*0(2) gauge transformation on the normal vectors in the Frenet frame 
(2.9): 

E 1 = E 1 = T, E 2 = E 2 cos6 + E 3 sin6, E 3 = -E 2 sin 6 + E 3 cos 9 (2.18) 
with the rotation angle 6 defined by 

O x = -r (2.19) 

so thus 

E lx = kcos6E 2 — k sin 9E 3 JL T, (2.20) 
E 2x = -KCOS6E! || T, E 3x = Ksm6E 1 \\ T. (2.21) 

This is called a parallel framing [7] of the space curve f(x). The frame vectors (2.18) are 
characterized by the geometrical property that along f{x) their derivatives lie completely in 
the normal space (2.20) or in the tangent space (2.21). Such a frame is unique up to a rigid 
(^-independent) rotation 

9 ->■ 6 + 0, = const. (2.22) 
acting on the pair of normal vectors. 
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In matrix notation the Serret-Frenet equations of a parallel frame are given by 



E X = UE 



(2.23) 



with 



/ 



E 



T 

cos ON + sin 6B 
- sin ON + cos 6B 



\ 



U 



V 



u 2 u 3 
—u 2 u\ 

— U 3 —Mi 



eso(3) 



(2.24) 



where 



Ui = 0, m 2 = kcosO = /tcos(J rrfx), u 3 = —KsinO = /tsin(J rdx) (2.25) 
are the components of the principal normal T x of f(x). The evolution of this frame 



E, = WE 



(2.26) 



is described by the matrix 



W = 



ZU 2 ZU3 
—ZU2 VJ\ 

-w% —w\ 

which is related to U through the zero-curvature equation 



V 



eso(3) 



(2.27) 



/ 



U t -W x + [U,W] = 0. 



(2.28) 



Note W can be determined directly from the model (2.1) via the relations (2.9) and (2.11). 
It now becomes convenient to work in terms of a complex variable formalism 



zu = ZU2 + IZU3, 

u = U2 + ius = ne~ ld = K,exp(i J rdx), 



(2.29) 
(2.30) 



encoding the well-known Hasimoto transformation [13]. Ex. in the Heisenberg model (2.2), 
the vortex filament equations on k and r transform into the NLS equation on u: 

1, 



- iu t = u xx + - \u\ u. 



(2.31) 



Thus, Hasimoto's transformation has the geometrical interpretation [10] of a SO (2) gauge 
transformation on the normal frame of the curve f(x), relating a Frenet frame to a parallel 
frame. 



Remark: Since the form (2.25) of a parallel frame is preserved by 5*0(2) rotations 
(2.22), the complex scalar variable (2.30) given by the Hasimoto transformation is uniquely 
determined by the curve f(x) up to rigid phase rotations u — > e~ l ^u, depending on an 
arbitrary constant <f). Therefore, u has the geometrical meaning of a covariant of the curve 
[5] relative to the group SO (2) ~ U(l), while \u\ = k and (&rgu) x = r are invariants of the 
curve. 

III. BI-HAMILTONIAN FLOWS AND OPERATORS 

For a general vector model (2.1) the zero-curvature equation (2.28) gives an evolution 
equation on u, 

u t — zu x — iw\U, (3.1) 
plus an auxiliary equation relating w\ to u, 

W\ x = lm(wu). (3.2) 

From (3.2) we can eliminate zu 1 = D~ 1 lm(zuu) in terms of u and zu, and then we see (3.1) 
yields 

u t = D x w — iuD~ l \m.{wu) = %{w) (3.3) 
where w is determined from (2.1) via the frame evolution equation (2.26). 

Proposition 1. 

U = D X - iuD- l lm(uC) (3.4) 

is a Hamiltonian operator with respect to the flow variable u(t,x), whence the evolution 
equation (3.3) has a Hamiltonian structure 

u t = H(5Sj/5u) (3.5) 

iff 

zu = 5H/5u (3.6) 

holds for some Hamiltonian 

fj = / H(x,u,u,u x ,u x ,u xx ,u xx , . . .)dx. (3.7) 
Jc 
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Here C is the complex conjugation operator, and C = K or S 1 is the domain of x. In the 
present setting, an operator V is Hamiltonian if it defines an associated Poisson bracket 

6} = / Re(V(8$j/8u)5<5/5u)dx (3.8) 

obeying skew-symmetry {$), (55} = — {<3, ij} and the Jacobi identity 0}} + cyclic = 0, 

for all real- valued functionals (3, on the x-jet space of the flow variable u. 

Proposition 2. (i) The Hamiltonian operator H is invariant with respect to U(l) phase 
rotations e iA "He~ iA = T-i\ u ~>e iX u- (H) A second Hamiltonian operator is given by 

T=-i (3.9) 

which is similarly U (1) -invariant, e lX Te~ %x = X. (Hi) The operators H and! are a compatible 
Hamiltonian pair (i.e. every linear combination is again a Hamiltonian operator), and their 
compositions define U(l) -invariant hereditary recursion operators 

Tl = Hl- 1 =i(D x + uD~ 1 Re(uC)), K* = T^U = iD x - uD- l Re(iuC). (3.10) 

(iv) Composition of 71 and H yields a third U '(1) -invariant Hamiltonian operator 

S = TZn = iD 2 x + D x (uDj 1 Im(uC)) + iuD- l Re{uD X C) (3.11) 
= iD 2 x + i\u\ 2 + u x D x l lm{uC) - iuD~ 1 'Re(u x C) 

satisfying e lX £e~ lX = £\ u ^ e i\ u . In particular, £, H, 1 form a compatible Hamiltonian triple. 

These Propositions are a special case of group-invariant bi-Hamiltonian operators derived 
from non-stretching curve flows in constant-curvature spaces and general symmetric spaces 
in recent work [3-5, 24]. Moreover, in the present complex variable formalism, Proposition 2 
provides a substantial simplification of some main results in [17] on Hamiltonian operators 
connected with non-stretching curve flows in Euclidean space. 

Because phase rotation on u is a symmetry of both % and X, the recursion operator 1Z 
generates a hierarchy of commuting Hamiltonian vector fields given by 

izu (n) d/du = 7Z n (iu)d/du, n = 0, 1, 2, . . . (3.12) 

where 

w (n) = 5HW/5u = K* n (u), n = 0,l,2,... (3.13) 
11 



are Hamiltonian derivatives, starting with 

w^=u, H^=uu=\u\ 2 (3.14) 
which corresponds to phase-rotation iud/du. Next in the hierarchy comes 

ujW = iu x , if (1) = -{uu x - uu x ) = lm(u x u), (3.15) 

followed by 

zu^ = -(u xx + ^\u\ 2 u), = \u x \ 2 - l -\u\\ (3.16) 

corresponding to respective Hamiltonian vector fields —u x d/du which is rc-translation and 
—i(u xx + \ \u\ 2 u)d/ du which is of NLS form. 

Through Propositions 1 and 2, this hierarchy produces integrable evolution equations 
on u(t,x) with a tri-Hamiltonian structure. An explicit formulation of this result has not 
appeared previously in the literature. 

Theorem 1. There is a hierarchy of integrable bi- Hamiltonian flows on u(t,x) given by 

u t = U(5rt n) /5u)=T(5rt n+1) /5u), n = 0,1,2,... (3.17) 

(called the +n How) in terms of Hamiltonians S)^ = j c H^dx where 

if W = -2—D x l \m{u{iU) n+1 u) n = 0, 1, 2, . . . (3.18) 

are local Hamiltonian densities. Moreover, all the flows for n ^ have a tri-Hamiltonian 
structure 

u t = £(5f) (n - 1) /5u), n = l,2,... . (3.19) 

Remarks: Each flow n = 0, +1, +2, . . . in the hierarchy is £/(l)-invariant under the 
phase rotation u — > e lX u and has scaling weight t — > \ 1+n t under the NLS scaling symmetry 
x — > \x, u — > \~ 1 u, where the scaling weight of ff (™) is —2 — n. Additionally, these flows on 
u(t,x) each admit constants of motion (under suitable boundary conditions) 

D t I \u\ 2 dx = 0, D t I iuu x dx = 0, D t f \u x \ 2 — -\u\ 4 dx = 0, ... (3.20) 
Jc Jc Jc 4 

and symmetries 

1 3 

-u x d/du, -i(u xx + -\u\ 2 u)d/du, (u xxx + -\u\ 2 u)d/du, ... (3.21) 
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respectively comprising all of the Hamiltonians (3.18) in the hierarchy and all of the corre- 
sponding Hamiltonian vector fields (3.12). 

At the bottom of the hierarchy, the flow is given by a linear traveling wave equation 
u t = u x , and next the +1 flow produces the NLS equation (2.31). The +2 flow yields the 
complex mKdV equation 

- u t = u xxx + ^\u\ 2 u (3.22) 
which corresponds to an mKdV analog of the vortex filament equations, 

— Kt — (k, xx + —k ) x — — = n xxx + — (ac — 2r )n x — 3ktt x , (3.23) 

Z Z K z 

-n = (r xx + 3 + \rK 2 - r% (3.24) 



\ T J^L + Tx (q1^ - 3 4 - 3r 2 + -K 2 ) + r(3^ - 3^^ 

K K K 2 K K 



as obtained through the Hasimoto transformation u = Kexp(-iO). 

The evolution equations describing the 0, +1, +2, . . . flows on u each arise from geomet- 
ric space curve flows corresponding to S'0(3)-invariant vector models (2.1). To make this 
correspondence explicit, it is convenient to introduce a complex frame notation 

E ll =E 1 = T, E 1 - = E 2 + iE 3 = e~ i0 (N + iB) (3.25) 

satisfying 

AE ± = -iE 1 - = E 3 - iE 2 = e~ i9 {B - iN) (3.26) 

and 

£ll.£ll = l, E ± -E ± = 2, E 11 ■ E L = = E L ■ E ± . (3.27) 
The Frenet equations (2.14) become 

£j = Re(^), E^ = -uE\ (3.28) 
while from (2.27), (2.29), (3.2), the evolution of the frame is given by the equations 

E\ = Re(wE ± ), = iD~ 1 lm(wu)E ± - wEK (3.29) 

Then any flow belonging to the general class 

zu = w(u, u, u x , u x , u xx , u xx , . . .) (3.30) 
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will determine a vector model (2.1) via 

S = E\\ f = Re(wE ± ), (3.31) 

where / is expressed in terms of S, S x , S xx , etc. through the Frenet equations (2.23)-(2.24). 
Ex. 1: The +1 flow w = iu x yields 

e\ = -Re(iu x E L ). (3.32) 

By rewriting 

u x E ± = {uE L ) x + uuE 11 

we obtain 

ReiiuE^ = -EW A Re(-ui?- L ) = -EW A El, 
Re{iuuEW) = Re(i\u\ 2 )EU = 0, 

and hence 

Ef = (E^AEl) x . (3.33) 

The identifications (3.31) then directly give the SO (3) Heisenberg model (2.2), which cor- 
responds to the non-stretching space curve flow (2.17) or equivalently 

f t = kB, \f x \ = l (3.34) 

expressed as a geometric flow. 

Ex. 2: The +2 flow w = —(u xx + \\u\ 2 u) yields 

- e\ = Re((u xx + ^ufu^) = Re(« sx £;- L ) + \\u\ 2 El (3.35) 

Here we can rewrite the first term as 

Re{u xx E L ) = Re{uE L ) xx - Re(uE^) x - Re(u x E^) = E\ xx + (\u\ 2 E% + ]^{uu) x E\ 

with \u\ 2 = \E X \ 2 , and thus 

-^^^L + ld^Il 2 ^ 11 ).- (3-36) 

Hence, — S gives 

3 2 

— S t = S xxx + — (l^xl S) x (3.37) 



14 



which can be viewed as an 5*0(3) mKdV model. The corresponding non-stretching space 
curve flow looks like 

3 

— Tf — r xxx + — \r xx \ r x , \r x \ = 1. (3.38) 
This describes a geometric flow [16] 

_f t = I K 2 T + k x N + ktB, \f x \ = l (3.39) 

which is equivalent to the evolution (3.23) and (3.24) on the curvature and torsion of f(x). 

Remark: A different geometric derivation of the mKdV model (3.37) appears in work [2] 
on non-stretching flows of curves in three-dimensional manifolds with constant curvature, 
i.e. S 3 , H 3 , R 3 , where the spin vector S is identified with the components of the unit 
tangent vector in a moving frame defined by parallel transport along the curve. The mKdV 
model also has been derived in [1 1] as a higher-order symmetry of the Heisenberg model by 
non-geometric methods. 

All of these SO (3) vector models describe dynamical maps 7 on the unit sphere S 2 C R 3 
by means of the identifications: 

S t ^ -ft, S x <->• 7x , d x + S{S X -) ^V x , SA^J = * (3.40) 

and thus 

22 2^2 

^ xlx S xx + I S x I S, V \fix ^ S xxx + \S X \ S x + — (\S X \ ) X S, (3-41) 
J lx ^S A S x , JV X -/ X ^ S A S xx , (3.42) 
g(lx,1x) = \lx\ 2 g ^ \S X \ 2 = S x - S x (3.43) 

where g denotes the Riemannian metric on the sphere S 2 (given by restricting the Euclidean 
inner product in M 3 to the tangent space of S 2 C M 3 ). 

In particular, the SO (3) Heisenberg model yields the Schrodinger map equation (2.4) on 
S 2 , while the SO (3) mKdV model (3.37) is identified with 

-lt = Vllx + ^\lx\)lx (3.44) 

which is a mKdV map equation on S 2 (i.e. a dynamical map version of the potential mKdV 
equation). 

Thus, Theorem 1 provides a geometric realization of the hierarchies of integrable vector 
models and dynamical maps containing the Heisenberg model and the Schrodinger map as 
well as their mKdV counterparts. 

15 



IV. GEOMETRIC HIERARCHY OF INTEGRABLE VECTOR MODELS AND 
DYNAMICAL MAPS 



In general, any non-stretching space curve flow (2.7) can be written in terms of a Frenet 
frame (2.9) by an equation of motion of the form 

f t = aT + bN + cB (4.1) 

such that 

D x a = Kb. (4.2) 

This relation between the tangential and normal components of the motion arises due to the 
non-stretching property 

\r x \ = 1 (4.3) 

by which the motion preserves the local arclength ds = \f x \dx of the space curve if (and 
only if) f x ■ f tx = 0. As a consequence, through the Serret-Frenet equations (2.14)-(2.15), 
the tangent vector T = r x along the space curve obeys the equation of motion 

T t = f\N + f 2 B = f JL T (4.4) 

given by a linear combination of the normal and bi-normal vectors with coefficients 

fx = D x b -tc + Ka, f 2 = D x c + rb. (4.5) 

Now we consider a Hasimoto transformation (2.18)-(2.19) from the Frenet frame (2.9) to 
a parallel frame (3.25). The equation of motion (4.1) on r takes the form 

f t = + Re(/i ±J B ± ) = h\\T + Re(h ± e id )N + lm(h ± e w )B (4.6) 

in terms of the tangential and normal components given by 

h ± = (b + ic)e~ w , h\\ = a, (4.7) 

with these components satisfying the relation (4.2) given by 

D x h\\ = Re(uh ± ) (4.8) 

where u = ne~ %e . Correspondingly, from the Frenet equations (3.28) of the parallel frame, 
the equation of motion (4.4) for the tangent vector T — r x has the form 

T t = Re(roE ± ) = Re(zue id )N + lm(zue ie )B (4.9) 
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in terms of 

zu = D x h± + h\\u (4.10) 
which encodes the normal and bi-normal components 

A + if 2 =zue ie . (4.11) 

The evolution of T is thus specified by the variable zu, while the underlying evolution 
of r is specified in terms of the variable h±, with h\\ given by the non-stretching condition 
(4.8). From equation (4.10) these variables are related by 

zu = D x h ± + uD^Reiuh^ = J(h ± ). (4.12) 

The operator here 

J = D x + uD^ReiuC) (4.13) 
is related to the Hamiltonian operator X = — i by the properties 

- J = TZ*!- 1 = X^TZ and - J' 1 = TZ~ X X = IK*' 1 , (4.14) 

where TZ and TZ* are the recursion operators (3.10). Consequently, J~ l is a formal Hamil- 
tonian operator compatible with X. 

Proposition 3. The evolution (4-1) of a non- stretching space curve f(x) can be expressed 
in terms of a geometrical variable that determines the corresponding evolution (4-4) °f the 
tangent vector T = r x through the relation T t = (f t ) x . In particular, 

h ± = J-\w) = TZ-\iw) = iTZ*-\w) (4.15) 

yields the normal components of the evolution vector f t in a parallel frame, where zu repre- 
sents the frame components of T t . The curvature k and torsion r of f(x) correspondingly 
have the evolution 

n t = D x f x - rf 2 = Re(e w D x zu) (4.16) 
r t = D x {K- x D x f 2 + tk- x K) + Kf 2 = D^K-'lmie^D^)) + dm.{e*w) (4.17) 

which can be expressed in terms of the Frenet frame coefficients a, b, c of f t through the 
relations (4-5). 
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Conditions will now be stated within the general class of flows (3.30) on u such that 
the various evolutions (4.16)-(4.17), (4.9), (4.6), (4.4), (4.1), and (3.5)-(3.6) each define a 
geometric flow. 

Theorem 2. For a non- stretching flow of a space curve f(x) in M 3 7 the following conditions 
are equivalent: 

(i) Its tangent vector T = f x = S obeys a SO (3) -invariant vector model iff fi and f 2 are 
functions of scalar invariants formed out of S and its x derivatives (modulo differential 
consequences of S • S = 1), i.e. 

/1+2/2 = f(S x -S x , S x -S xx , S xx -S xx , . . . , S xx -{SaS x ), S xxx -(SaS x ), S xxx -(SaS xx ), . . .). (4.18) 

(ii) Its principal normal component u = T x ■ E 1 - in a parallel frame (3.25) satisfies a U{1)- 
invariant evolution equation iff w is an equivariant function of u, u, and x derivatives of 
u and u, under the action of a rigid (x -independent) U(l) rotation group u — > e~ l<i> u (with 
4> =const.), i.e. 

™ = uf{\u\, \u\ x , \u\ xx , . . . , (argw)^ (argu) xx , . . .). (4.19) 

(Hi) Its curvature k and torsion r satisfy geometric evolution equations in terms of invariants 
and differential invariants off(x) iff we ld is a function of k, t, and their x derivatives, i.e. 

ZU = e K-f^K, K X) K XX) . . . , T, T X) T xx , . . .). (4.20) 

(iv) Its equation of motion is invariant under the Euclidean isometry group SO(3) x R 3 iff 
a, b, c are scalar functions of the curvature k, torsion r, and their x-derivatives, subject to 
the non- stretching condition (\.2). 

The proof of this proposition amounts to enumerating the Euclidean (differential) invari- 
ants of a space curve r[x) with an arclength parameterization x, as shown in appendix A. 

We are now able to derive the entire hierarchy of S'0(3)-invariant vector models and 
geometric space curve motions that correspond to all of the [/(l)-invariant flows on u in 
Theorem 1. 

From equation (4.9) combined with the hierarchy (3.6), the evolution of the spin vector 
S — T — f x can be written as 

S t = Re(E ± IZ* n (u)), n = 0,l,2,... (4.21) 
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as generated via the recursion operator TV = iD x — uD x 1 Re(iwC). The main step is now to 
establish the operator identity 

Re(E L TZ*) = SRe^C) (4.22) 

where S is a spin vector operator corresponding to TV, and Re(E L C) is the operator that 
produces a vector / in the perp space of S in IR 3 when applied to the components of / with 
respect to E L (i.e. / = Re{E L Cf) if / • S = 0, where / = / • E L ). Through the Frenet 
equations (3.28) and the orthonormality relations (3.27) on E 1 - and E^, we straightforwardly 
find 

Re(E ± u) = El = S x (4.23) 

and 

Re(E L iDJ) = A Re(E L D x f) = S A D x Re(E L Cf) = S A D x f 
D-'ReiiuCf) = D^Re^ A E\) ■ E^Cf) = D^Re^S A S x ) ■ f) 

for any vector f(x), orthogonal to S in IR 3 , with components f(x) = E L • f(x). Hence, this 
yields the vector operator 

S = S A D x - S X D7 £ \S A S x )- (4.24) 

Theorem 3. The bi-Hamiltonian flows (3.17) on u(t,x) correspond to a hierarchy of inte- 
grate SO (3) -invariant vector models 

S t =(SAD x -S x D-\SAS x y) n S x = f (n \ n = 0,1,2,... (4.25) 

generated by the recursion operator (4-24)- These models have the equivalent geometrical 
formulation 

lt={JV x - lx D^g{J lx , )) n lx = F^, n = 0,l,2,... (4.26) 

expressed as dynamical maps j(t,x) into the 2-sphere S 2 CK 3 . Moreover, the Hamiltonians 
(3.18) for all the flows on u(t, x) correspond to a set of constants of motion = J c H^dx, 
ijW = J c H^dx, S)^ = J c H^dx, etc. for each vector model (4-25) and each dynamical 
map equation (4-26). This entire set has the explicit form given by the densities (modulo 
total x- derivatives) 

(1 + n)HW = D7 X \S X • D x /W) = D-'g^, V X F™), n = 0, 1, 2, . . . (4.27) 
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which are scalar polynomials formed out of SO (3) -invariant wedge products and dot products 
of S, S x , S xx , . . . in terms of the equations of motion for S(t, x), or equivalently, scalar inner 
products of~f x , J^ x , V x j x , JVxIx, ■ ■ ■ formed in terms of the equations of motion for j(t, x). 

The vector models (4.25) have been derived previously in [6] by non-geometric methods 
(based on a Lax pair representation). As one new result, Theorem 3 derives the equiva- 
lent dynamical map equations (4.26), along with their recursion operator, and provides an 
explicit expression for the constants of motion for all of these integrable generalizations of 
the Heisenberg spin model (2.2) and the mKdV spin model (3.37). In particular, from the 
hierarchy (4.25), higher-derivative versions of the Heisenberg spin model n = 1 as given 
by n — 3, 5, . . . are seen to describe higher- derivative Schrodinger maps; similarly, higher- 
derivative versions of the mKdV spin model n = 2 as given by n — 4, 6, . . . are found to 
describe higher- derivative mKdV maps. 

Ex: n = 3 and n = 4 respectively yield a 4th order Heisenberg vector model 

S t = S A S x X xx + \{\Sx\ 2 S A S x )x (4.28) 
and a 5th order mKdV vector model 

— St = Sxxxxx — -^\S XX \ S + \S X \ S x — — \S X \ S) x + — (\S X \ S)xxx (4.29) 
which are described geometrically by a 4th order Schrodinger map equation 

- 7* = J^llx + ^(17*1^7*) - \9{Jlx, V xlx )lx (4.30) 
and a 5th order mKdV map equation 

It = ^tlx + ^{\lx\ 2 g)xV x lx + (^(iTxIfl)^ - ^Xlx\l + ^\lx\g)lx- (4.31) 

Remark: The equations of motion of these dynamical maps 7(t, x) do not locally preserve 
the arclength ds = \^ x \gdx in the x direction along 7 on S 2 , namely |7 x | g has a nontrivial 
time evolution for any of the dynamical map equations (4.26), and additionally, the total 
arclength f c \^ x \ g dx is time dependent. 

The spin vector recursion operator (4.24) has the factorization 

S = (D x + SS X ■ +S x D- 1 Sx- )(SA ) (4.32) 
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where, as shown by the results in [6], the operators D X +SS X - +S X D~ 1 S X - and (SA ) _1 = — SA 
constitute a compatible Hamiltonian pair with respect to the spin vector variable S(t, x). By 
comparison, using the present geometric framework, we now directly derive the explicit bi- 
Hamiltonian structure of the hierarchy of vector models (4.25) through the bi-Hamiltonian 
flow equation (3.17) on u(t,x) in Theorem 1. 
The starting point is the variational identity 



holding modulo total x-derivatives and modulo (differential consequences of) S ■ 8S — 0, 
with the Hamiltonian densities given in the equivalent forms (3.18) and (4.27), and with 
w ( n ) given by the hierarchy (3.13) generated through the recursion operator 1Z*. To begin 
we derive an explicit expression for 5u in terms of 5S = 5E^ by taking the Frechet derivative 
of the Frenet equation (3.28) for the normal vectors in a parallel frame with respect to u: 



The E 1 - component of (4.34) gives 

D^E 1 - ■ 5E L ) = -uE 1 - ■ SEW - uE^ ■ 5E ± = -ZiR^iuE 1 - • 5E 11 ) 

via E\\ ■ SE 1 - = —E 1 - ■ 5E^ due to the orthogonality of E L , EK Similarly, the component 
of (4.34) yields 



5H^ = 5S ■ (8HW/8S) = 2Re(w^5u) 



(4.33) 



D X SE ± = -5uE 11 -u5EK 



(4.34) 




Combining these two expressions, we obtain 



5u = D X (E 



■ 5E 11 ) - iuD^ReiiuE 1 - ■ SE 11 ) 



(4.35) 



whence 



Re(zu in) 5u) = Re(tu^D x (E ± ■ <?£")) + Im(ro (n) M)D; 1 Re(wE ± • SE^). 



(4.36) 



Integration by parts on both terms in (4.36) then yields 



Re(w in) Su) 



Re(E ± ■ 5EWD x w {n) ) - Re(iuE L ■ 5E^D~ 1 lm(uj in) u)) 
SEW ■ Re(E ± ^(w (n) )) 



(4.37) 
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modulo total x-derivatives. Next, after use of the relations ilZ* = —H and iE 1 - — A E 1 
which follow from (3.10) and (3.26), we get 

Re{w {n) 5u) = 60 ■ (£" A Re^ro^ 1 ))). (4.38) 

Hence, equating (4.38) and (4.33) yields 

SS ■ (5H {n) /5S - 2S A Re(E ± a7 (ri+1) )) = 

from which we obtain the variational relation 

- S A (5H (n) /5S) = 2Re(E ± ^ n+1) ) = 2Re(E ± n* n+1 (u)) = 2f {n+1) (4.39) 

where the final equality comes from the equation of motion (4.21) combined with the hier- 
archy (4.25). This result (4.39), together with the factorization of the recursion operator 
(4.32), leads to the following Hamiltonian structure. 

Theorem 4. In terms of the Hamiltonian densities (4-27), the hierarchy of vector models 
(4-25) for S(t,x) has two Hamiltonian structures 

S t = -S A (SH^-V/SS) = f in \ n = l,2,... (4.40) 

and, for n^l, 

S t = D x {5H {n ~ 2) /5S + SD-\S X ■ 5H (n ' 2) /5S)) = f {n \ n = 2, 3, . . . (4.41) 

where 

S A and D X + SS X - +S x D~ l S x - (4.42) 

are a compatible pair of Hamiltonian operators. Correspondingly, the hierarchy of dynamical 
maps (4-26) on j(t,x) has the Hamiltonian structures 

lt = -J(5H {n ~ 1) /5 1 ) = F {n \ n = l,2,... (4.43) 
= V^SH^/S-y) + i x D~ x g{ lx , SH^/S-y), n = 2, 3, . . . (4.44) 

in terms of the compatible Hamiltonian operators 

J and V x + lx D^ l g( lx , ) (4.45) 

as given by the geometrical identifications (3.40). Moreover, these bi- Hamiltonian pairs 
(4-42) and (4-45) are geometrically equivalent to the pair of compatible symplectic (inverse 
Hamiltonian) operators — X -1 and J = —TZ*!^ 1 with respect to the flow variable u(t,x) in 
the hierarchy (3.17)-(3.18) (cf. Theorems 1 and 3). 
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These results provide an explicit geometrical formulation of the abstract symplectic struc- 
ture given in [28] for the Schrodinger map equation (2.4), i.e. (n = 1) 

It = JV xlx = -J(6HW/6>y), # (0) = ^7(7*,7*), (4-46) 

and its higher-order generalizations (n = 3, 5, . . .); the Hamiltonian structure of the mKdV 
map, i.e. (n = 2) 

It = ~V 2 xlx - \\ lx \ 2 glx = V^SH^/S-y) + TxD-'gir,, 6H^/&y) (4.47) 
= -J(8HM/8 1 X H^ = ~g{V xlx ,J lx ), (4.48) 

and its higher-order generalizations (n — 4, 6, . . .) has not appeared previously in the litera- 
ture. Note the n = 1 and n = 2 vector models respectively yield the well-known Hamiltonian 
structure of the Heisenberg spin model (2.2), 

S t = S A S xx = -S A (5H^/5S), # (0) = ^\S X \ 2 , (4.49) 

and the mKdV spin model (3.37), 

S t = S xxx - ^(\S X \ 2 S) X = D X (6H<V/6S + SD^(S X ■ SH^/SS)) (4.50) 

= —S A (6HW/6S), H (1) = l -S x ■ (S A S xx ). (4.51) 

Remark: There is a second Hamiltonian structure for both the Schrodinger map equation 
and the Heisenberg spin model. This structure, in contrast to the first Hamiltonian structure 
(4.46) and (4.49), turns out to involve a non-polynomial Hamiltonian density defined as 
follows. Let £(7) be a vector field on S 2 such that its divergence is constant at all points 
7 G S 2 . (This property geometrically characterizes £(7) as a homothetic vector with respect 
to the metric-normalized volume form e g on S 2 , i.e. C^e g = ce g for some constant c ^ 0.) 
Then the Hamiltonian density given by 

H ( - l) = <?(£( 7 ), JV 7 «), div s £(7) = 1 (4.52) 

can be shown to satisfy (see Appendix B) 

5H^ =g(5 1 ,JV lx ) (4.53) 

modulo total x-derivatives, so thus 

lt = V x (6H(-V/6>y) + i x D~ l g( lxi SH^/S-y) = V x (J lx ) (4.54) 
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yields a second Hamiltonian structure for the Schrodinger map (2.4). The corresponding 
second Hamiltonian structure for the Heisenberg spin model (2.2) is given by 

S t = D X (6H(-V/6S + SD-\S X ■ SH^/SS)) — (S A S x ) x , — £(S) ■ (S A S x ) (4.55) 

in terms of a vector function £(S) such that 

S-t(S) = 0, d£-t(S) = l, (4.56) 

where the operator d$ = ds — S(S ■ ds) is the orthogonal projection of a gradient with 
respect to the components of S. From the properties 

S-dj = 0, dj{S-S) = 0, (4.57) 

the Hamiltonian density can be shown to satisfy 

5H { ~ 1] = 5S ■ (S A S x ) (4.58) 

modulo total x- derivatives (see Appendix C). Thus, both the Heisenberg spin model and 
the Schrodinger map equation are Hamiltonian equations of motion with respect to the 
corresponding bi-Hamiltonian pairs (4.42) and (4.45). 

To conclude, a counterpart of Theorem 3 will now be stated for the underlying space 
curve motions on r. 

Theorem 5. The SO (3) -invariant vector models (4-25) on S(t,x) correspond to a hierarchy 
of integrable flows of non- stretching space curves 

f t = a {n - 1) T + b {n - 1) N + c {n - 1) B, \f x \ = l, n = l,2,... (4.59) 

with geometric coefficients 

= Re(Q n K), = -Im(g n re), (4.60) 

a W = -D- 1 (/cIm(Q n K)) (4.61) 

given in terms of the recursion operator 

Q = iD x — t — kD-^kIui) = e ie TZ*e- ie , (4.62) 

where k,t are the curvature and torsion off. The tangential coefficients (4-61) yield a set 
of non-trivial constants of motion J c a^dx = — J c ^K 2 dx, J c a^dx = j c TK 2 dx, etc. for 
each curve flow (4-59) in the hierarchy (under suitable boundary conditions) , where C = R 
or S 1 is the coordinate domain of x. 
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The equations of motion (4.59) arise from writing the flow equation (4.6) on r in terms 
of the hierarchy (3.6) by means of the relation (4.15). This yields 

a {n) = L>; 1 Re(mw (n - 1) ) = -D^lmiuw^), (4.63) 

b {n) = Re(ie t6 w {n ~ l) ) = -lm(e u> w ( - n - 1) ), (4.64) 

c (n) = \ui(ie ie w {n ~ l) ) = Re(e i0 w {n - l) ). (4.65) 

The geometric form (4.60)-(4.61) for these coefficients is then obtained through the identity 

e ie w {n-i) = e iO n *n-i u = ^ e n*e' ie ) n - l e ld u (4.66) 

combined with the Hasimoto transformation 

K = e i9 u, t = -9 X . (4.67) 



V. SURFACES AND SPATIAL HAMILTONIAN CURVE FLOWS 

The results in Theorem 1 for 1+1 flows can be generalized in a natural way to 2+1 flows 
by considering surfaces f(x, y) that are foliated by space curves with a parallel framing in 
IR 3 . Here y will denote a coordinate assumed to be transverse to these curves, and x will 
denote the arclength coordinate along the curves, so thus 

\f x \ = 1. (5.1) 

In this setting, a parallel frame consists of a triple of unit vectors whose derivatives along 
each y = const, coordinate line on the surface f(x,y) lie completely in either the tangent 
space or the normal plane of this line in R 3 . The explicit form of such a frame is given by 
the vectors 

E 1 = T = f x , (5.2) 
E 2 = Re{e- W {N + iB)) = ^^(costf f xx + sin# f x A f xx ), (5.3) 
E 3 = lm(e- ie (N + iB)) = |f M | _1 (-sin0 f xx + cos9 r x Ar xx ), (5.4) 



with 



@x |^xx| (Tx ^ T'xx) ' T'xxx (5-5) 
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where T, N, B respectively denote the unit tangent vector, unit normal and bi-normal vectors 
of the y = const, coordinate lines. 

In matrix notation these frame vectors satisfy the Frenet equations 



E x = UE, E y = VE 



(5.6) 



given by 



E = 



where 



E 2 
\E 3 j 



( u 2 u 3 ^ 



, U 



-u 2 
-« 3 



^ v 2 v 3 ^ 



G50(3), V = 



7 



—v 2 vi 
-v 3 -Vi 



eso(3), (5.7) 



u 2 + iu 3 = u = (E 2 + iE 3 ) ■ E lx , v 2 + iv 3 = v = (E 2 + iE 3 ) ■ E ly (5.8) 
are complex scalar variables, and 

Vl = -(E 3 + iE 2 ) ■ (E 2y + iE 3y ) = ^ (5.9) 

is a real scalar variable. Note, similarly to the case of space curves, here we have 

(E 3 + iE 2 ) ■ (E 2x + iE 3x ) = (5.10) 

which characterizes E(x, y) as a parallel frame adapted to the foliation of the surface f(x, y) 
by re-coordinate lines, i.e. E\ x _L T, E 2x || T, E 3x || T. From (5.5) we see this choice 
of framing is geometrically unique up to rigid rotations that act on the normal vectors 
(5.3)-(5.4): 



9 ->■ 9 + 0, E 2 + ?-E 3 -> exp(-i0)(£ , 2 + i-E 3 ), with =const. . 
Under such rotations, v\ is invariant, while u and f undergo a rigid phase rotation, 



(5.11) 



(5.12) 



To proceed we need to write down the tangent vectors of f(x, y) in terms of the parallel 
framing, 

r x = E u f y = q 1 E 1 + q 2 E 2 + q 3 E 3 . (5.13) 
The respective projections of r y orthogonal and parallel to r x are given by the scalar variables 



Q2 + iQ3 = q = (E 2 + iE 3 ) ■ f y , q 1 =E 1 - f y , 



(5.14) 
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where, under rigid rotations (5.11) on the normal vectors in the parallel frame, q\ is invariant 
while q transforms by a rigid phase rotation, 

q ->■ e-^q. (5.15) 

Remark: Due to the transformation properties (5.12) and (5.15), the variables u, v, q 
represent U {l)-covariants of f(x, y) as geometrically defined with respect to the x-coordinate 
lines, where U(l) is the equivalence group of rigid rotations (5.11) that preserves the form 
of the framing (5.2)-(5.5) for the surface f(x,y). These variables will be seen later (cf. 
Propositions 7 and 8) to encode both the intrinsic and extrinsic surface geometry of f(x, y) 
in W 3 . 

We will now show that all surfaces f(x, y) with the x coordinate satisfying the non- 
stretching property (5.1) have a natural geometrical interpretation as a spatial Hamiltonian 
curve flow with respect the y coordinate. This interpretation arises directly from the struc- 
ture equations satisfied by the parallel frame (5.2)-(5.4) and the tangent vectors (5.13) 
adapted to these coordinates. Firstly, the frame connection matrices given in (5.7) obey the 
zero-curvature condition XJ y — V x + [U, V] = 0, which yields the structure equations 

u y — v x + iv\u = 0, v\ x — lm(vu) = 0. (5.16) 

Secondly, the frame expansions of the tangent vectors given by (5.13) obey the zero-torsion 
condition (f x ) y = (r y ) x , leading to the structure equations 

v-q x ~qiu = 0, q lx - Re(qu) = 0. (5.17) 

Proposition 4. Through the frame structure equations (5.16) and (5.17), the U(l)- 
covariants u, v, q of r(x, y) are related by the Hamiltonian equations 

u y = n(v), v = J(q), (5.18) 

with J = —X^IZ = —TZ*!^ 1 , where T-L = D x — iuD'Hm^C) and I — — i are a pair of 
compatible Hamiltonian operators with respect to u, and 1Z = "HZ -1 = i(D x + uD^Re^uC)) 
and 1Z* = I~ lr H = iD x — uD^Re^uC) are the associated recursion operators. 

This leads to a main preliminary geometric result. 
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Lemma 1. Let u be an arbitrary complex-valued function of x andy, and define 

v = n~ 1 (u y ) = in' 1 ( Uy ), q = j-\v) = -n- 2 ( Uy ), (5.19) 

Vl = D^lmivu), gi = Re(qu) (5.20) 

in terms of the function u and the formal inverse operators 7-L~ x and J~ l , where TZ~ l is the 
formal inverse of the recursion operator 1Z — Hi — ij. Then the matrix equations (5.6) 
and the vector equations (5.13) constitute a consistent linear system of 1st order PDEs that 
determine a surface f(x,y) with a parallel framing of the y = const, coordinate lines for 
which x is the arclength. 

For a given function u(x,y), the resulting surface f(x,y) and frame E(x,y) are unique 
up to Euclidean isometries. 

VI. SURFACE FLOWS AND 2+1 SOLITON EQUATIONS 

The derivation of bi-Hamiltonian soliton equations from non-stretching space curve flows 
in Theorem 1 and their geometrical correspondence to integrable vector models and dynam- 
ical map equations in Theorems 3 and 4 will now be generalized to geometrically analogous 
surface flows in M 3 . This will mean we consider surfaces that are non-stretching along one 
coordinate direction yet stretching in all transverse directions. (No assumptions will be place 
on the topology of the surface.) 

Such surface flows f(t, x, y) can be written naturally in terms of a parallel frame (5.2)-(5.4) 
adapted to the non-stretching coordinate lines, for which x will be the arclength coordinate 
and y will be a transverse coordinate, by an equation of motion 

r t = h 1 E 1 + h 2 E 2 + h 3 E 3 , \f x \ = l, (6.1) 

with 

h = D- l Re{uh), h 2 + ih 3 = h. (6.2) 

Here hi and h respectively determine the components of the flow that are tangential and 
orthogonal to the non-stretching x-coordinate lines on the surface. The relation (6.2) imposes 
the non-stretching property by which the flow (6.1) preserves the local arclength ds = \r x \dx 
along these coordinate lines, with u = u 2 + iu 3 given by the components of the parallel frame 
connection matrix with respect to x from (5.7). 
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The corresponding evolution of the parallel frame is given by the same matrix equations 

(2.26) -(2.27) that govern a parallel framing of non-stretching flows of space curves. Equiv- 
alently, in the complex variable notation (2.29) for the components of the evolution matrix 

(2.27) , the evolution equations on the frame vectors consist of 

E u = Re(w(E 2 - iE 3 )), E 2t + iE 3t = -wE x + w x {E 3 - iE 2 ), (6.3) 

where 

w = D x h + Wiu, w lx = Re(uh) (6.4) 
are given in terms of h and u. These equations have the following Hamiltonian interpretation. 

Proposition 5. The respective evolution equations (6.1) and (6.3) for the surface and the 
parallel frame are related by 

h = J~\w) =TZ- 1 (iw) =iTZ*-\w) (6.5) 

where J~ l is a formal Hamiltonian operator compatible with the Hamiltonian pair H and 
T, and where Tlr 1 , IZ*^ 1 are inverse recursion operators (cf. Proposition A). 

Hence, surface flows of the type (6.1)-(6.2) can be expressed in terms of the complex scalar 
variable w. This variable also determines the resulting evolution of the frame connection 
matrices (5.7) through the pair of zero-curvature equations 

U t - W x + [U,W] =0 (6.6) 
V t -W y + [V,W]=0 (6.7) 

which express the compatibility between the Frenet equations (5.6) and the evolution equa- 
tions of the parallel frame in matrix form (2.26)-(2.27). From (6.6) we see u satisfies the 
same evolution equation (3.3) in terms of w as holds for non-stretching space curve flows in 
Proposition 1. We then find that the evolution equation obtained from (6.7) holds identi- 
cally as a consequence of Lemma 1. This result leads to the following Hamiltonian structure 
derived from surface flows (6.1)-(6.2). 

Lemma 2. The Hamiltonian structure (3.5)-(3.7) for 1+1 flows generalizes to 2+1 flows 
on u(t, x, y) given by Hamiltonians of the form 

/ / H(X, y, M, It, U X , U X , Uy, Uy, U XX , U XX , U X y, U X y, Uyy, Uyy, . . .) AX 0\j (6.8) 

J Jc 
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(where C denotes the coordinate domain of(x,y)). In particular, 

zu = 5H/5u (6.9) 
yields a 2+1 Hamiltonian evolution equation 

u t = U{w), U = D X - iuD-HmiuC), (6.10) 
corresponding to a surface flow given by (6.1), (6.2), (6.5). 

Since the Hamiltonian operator % does not contain the y coordinate, it obviously has y- 
translation symmetry. Hence starting from —u y d/du, there will be a hierarchy of commuting 
Hamiltonian vector fields 

iw [n) d/du = -n n (u y )d/du, n = 0,1,2,... (6.11) 

where 

w (n) =5H (n) /5u = n* n (iu y ), n = 0,1,2,... (6.12) 

are derivatives of Hamiltonian densities . An explicit expression for these densities can 
be derived by applying the scaling symmetry methods in [1, 5]. 

Theorem 6. The recursion operator 71 = "HZ -1 = i(D x +uD~ 1 Re(uC)) produces a hierarchy 
of integrable bi- Hamiltonian 2+1 flows 

u t = H{5& n - 1) /u) =l{5f) {n) /5u), n = l,2,... (6.13) 

(called the +n flow] in terms of the compatible Hamiltonian operators H and X, with the 
Hamiltonians S)^ = JJ C H^dxdy given by 

H (n) = Y^D-^uK^tuy)), n = 0,l,2,... (6.14) 

(modulo total x,y -derivatives). 

At the bottom of this hierarchy, 

H {0) = Re(iuu y ), 5H i0) /5u = iu y = ro (0) (6.15) 

yields the +1 flow 

- iu t = u xy + u u, u 0x = \u\\u\ y . (6.16) 
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This is a 2+1 nonlocal bi-Hamiltonian NLS equation which was first derived from Lax pair 
methods by Zhakarov [29] and Strachan [26]. 
Next in the hierarchy is 

H {1) =Re(u x u y ) -\uq\u\ 2 , 5H {1) /5u = -{u xy + u u) = w {1) . (6.17) 

This yields the +2 flow 

- u t = u xxy + \u\ 2 u y + u u x + iuiu, u lx = lm(u y u x ) (6.18) 

which is a 2+1 nonlocal bi-Hamiltonian mKdV equation. It can be written in the equivalent 
form 

-u t = u xxy + (u u) x + iv 2 u, v 2x = lm(uu xy ) (6.19) 

studied in work of Calogero [8] and Strachan [27]. 

These 2+1 flow equations have the following integrability properties. 

Proposition 6. The hierarchy (6.13)-(6.14) displays U{l)-invariance under phase rotations 
u — > e tX u and homogeneity under scalings x — > Xx, y — > Xy, u — > X~ 1 u, with t — > X 2+n t for 
the +n flow, where the scaling weight of is —3 — n. Each of the evolution equations 
(6.13) in this hierarchy admits the constants of motion 

D t J J iuu y dxdy = 0, D t J J u x u y — -vo\u\ 2 dxdy = 0, ... (6.20) 

and 

D t J J \u\ 2 dxdy = 0, D t J J iuu x dxdy = 0, D t J J \u x \ 2 — -\u\ 4 dxdy = 0, ... (6.21) 

comprising, respectively, all of the 2+1 Hamiltonians (6. 14) P^ us the 1+1 Hamiltonians 
(3.18) extended to two spatial dimensions (under suitable boundary conditions depending 
on the coordinate domain C of (x,y)). Additionally, these evolution equations (6.13) each 
admit the corresponding Hamiltonian symmetries 

— Uyd/du, —i{u xy + i>Qu)d/du, ... (6.22) 

plus 

iud/du, —u x d/du, — i(u xx + -\u\ 2 u)d/du, .... (6.23) 
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We note the constants of motion (6.21) and symmetries (6.23) are inherited from the 
1+1 integrability properties in Theorem 1 as a consequence of the fact that the Hamilto- 
nian phase-rotation vector field iud/du (which generates the hierarchy of 1+1 flows (3.17)) 
commutes with the Hamiltonian y-translation vector field —u y d/du (which generates the 
hierarchy of 2+1 flows (6.13)). 

VII. 2+1 VECTOR MODELS AND DYNAMICAL MAPS 

Each evolution equation (6.13) in the hierarchy presented in Theorem 6 determines a 
surface flow f(t, x, y) and a corresponding 2+1 vector model for S(t, x, y) = r x = E\ through 
the frame evolution equations (6.3) in a similar manner to the derivation of flows of space 
curves and 1+1 vector models. 

The +1 flow yields the geometric 5*0(3) vector model known as the M-I equation [19, 20, 

22] 

S t = SAS xy + v 1 S x = (SASy + v 1 S) x , v\ x = —S ■ (S x A S y ), (7.1) 

which is a 2+1 integrable generalization of the SO (3) Heisenberg model. It corresponds to 
the surface flow 

f t = r x A f xy + v x r x , v lx = -f x ■ (f xx A f xy ), \f x \ = I. (7.2) 

This flow equation describes the motion of a sheet of non-stretching filaments in Euclidean 
space, in analogy with the form of the vortex filament equations (2.8). Some properties of 
the model (7.1) have been studied recently in [9, 30]. 

The +2 flow produces a 2+1 integrable generalization of the geometric SO (3) mKdV 
model, 

— St — S xxy + ((S x • S y + vo)S — viS A S x ) x , vox = \S x \\S x \ y , (7.3) 
(the so-called M-XXIX equation [22]) which describes the surface flow 

T% T' XX y + (t xx ■ T' X y + VQ^jT' x V\T x A T xx , Vq x | T xx 1 1 T xx | y , \fx\ 1- C^ - ^) 

Each of these surface flows r(t, x, y) in IR 3 geometrically corresponds to a dynamical map 
7(i, x,y) on the unit sphere S 2 C M 3 through extending the identifications (3.41) and (3.42) 
as follows: 

Sy <-> 7y , 8y + S(Sy-) ^ Vy, (7.5) 
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S xy + (S x ■ S y )S <-> V xly = V ylx , (7.6) 

S A S xy +> JV x 7j, = JVyJ x . (7.7) 

The +1 and +2 flows thereby yield, respectively, 

It = J^xly + Ul7x, ui* = ^(7x, ^7y)> (7-8) 

and 

-7* = V^ + ^l^-ViJV^-^*, v 2x = g(i y , V^), (7.9) 

which are new nonlocal 2+1 integrable generalizations of the Schrodinger map equation (2.4) 
on S 2 and the mKdV map equation (3.44) on S 2 . 

The complete hierarchy of vector models and dynamical map equations in 2+1 dimensions 
can be written down in the same manner as in 1+1 dimensions (cf. Theorems 3 and 4) by 
means of the spin vector recursion operator (4.24) and its Hamiltonian factorization (4.32). 
In particular, the obvious y-translation invariance of this operator provides the geometric 
origin for the 2+1 generalization of the Heisenberg spin model and the mKdV spin model. 

Theorem 7. (i) The hi- Hamiltonian flows (6.13) correspond to a hierarchy of integrable 
SO '(3) -invariant vector models in 2+1 dimensions 

S t = (SAD x - S X D-\S A S x y ) n Sy = /W n = 1, 2, . . . (7.10) 

which are geometrically equivalent to 2+1 dimensional dynamical maps 7 into the 2- sphere 
S 2 c R 3 

lt={JV x - lx D- 1 g{J lx ,)) n ly = F^\ n = l,2,.... (7.11) 

(ii) Each vector model and dynamical map in the hierarchy (7.10)-(7.11) possesses a set 
of polynomial constants of motion that correspond to all of the Hamiltonians (6.14) f or the 
+l,+2,... flows (6.13), i.e. ^ = ff c H<®dxdy, = ff c H^dxdy, etc., as obtained 
from the Hamiltonian densities (modulo total x,y- derivatives) 

(1 + n)HM = -D-\S X ■ D x fW) = -D^g( lx , V X F^), n = 0, 1, 2, . . . (7.12) 

given in terms of the equations of motion (7.10) for S(t,x) and (7.11) for j(t,x). In addi- 
tion, the vector models (7.10) and dynamical maps (7.11) each possess two non-polynomial 
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constants of motion & ^ = ff c H^ ^dxdy and & 2 ) = ff c H( 2 ^dxdy explicitly given by 

H ( - 2) = \^S) ■ (S A S y ) = ^(£(7), Jly), (7.13) 
H [ - l) = \s x ■ S y + ^(S) ■ (S A S x ) = \g{ lx , ly ) + \v ig {t{i), J-Yx), (7-14) 

where £(7) is a vector field with covariantly- constant divergence div 9 £(7) = 1 at all points 
7 G S 2 7 and where ^(S) is an analogous vector function satisfying S ■ £(S) = 0, d§ ■ £(S) = 
1, in terms of the component-wise gradient operator 83 = 83 — S(S ■ 83) with properties 
(4-57). These Hamiltonian densities (7.13)-(7.14) correspond to two compatible nonlocal 
Hamiltonian structures for the +0 flow 

u t = £(5&-V/u) = n(8&-V/8u) = -u y (7.15) 

with ( cf. Lemma 1 ) 

- SH^-V/u = TZ^iiUy) = v, -5H ( - 2) /u = TZ*- 2 (iu y ) = -iq. (7.16) 

(Hi) In terms of the Hamiltonian densities (7.12), (7.13) and (7.14), all the 2+1 vector 
models (7.10) and dynamical map equations (7.11) have the bi- Hamiltonian structure 

S t = -S A (5H (n - 2) /5S) = D x {5H (n - 3) /5S + SD;\S X ■ 5H^/5S)) (7.17) 
= /<»>, n = l,2,... 

and 

lt = - J^H^/S-y) = V x (SH^/5j) + i x D- x g{ lx , SH^/S-y) (7.18) 
= F {n \ n = 1,2,... 

given by the respective pairs (4-4%) an d (4-45) of compatible Hamiltonian operators. 

Remark: Explicit bi-Hamiltonian formulations for the 2+1 generalization of the Heisen- 
berg spin model (7.1) and the geometrically corresponding new 2+1 integrable Schrodinger 
map (7.8) are given by 

St — —S A (-S xy + Vl S A S x ) = D X (S A S y + SD-\S X ■ (S A S y ))) (7.19) 

and 

7t = -J{-V y l x + «i J7x) = V x ( J 7l/ ) + IxD^gbx, Jl v ) (7.20) 
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where 

5H^ = 5S-(SA S y ) = g{5 1 , J ly ) (7.21) 
6H(-V = 5S ■ (S xy + Vl S A S x ) = g(Sj, -V y7 x + v 1 J lx ) (7.22) 

yield the respective derivatives of the non-polynomial Hamiltonian densities (7.13) and 
(7.14). These two densities in addition to all the polynomial densities (7.12) give a set 
of constants of motion for the Hamiltonian equations (7.19) and (7.20). In particular, the 
first four constants of motion are explicitly given by the integrals 

A JJj(S)-(SAS v )dxdy = D t JJ g^{i),Ji v )dxdy = Q, (7.23) 
A J J S x -S y + v^(S) ■ (S A S x )dxdy 

= A J J girix, ly) + v ig (Z(~f), Jl x )dxdy = 0, (7.24) 
A ^ • (S x A S^) + vi\S x \ 2 dxdy 

= D tJJ g{^ylx, Jlx) + vMldxdy = 0, (7.25) 
A // c • S„ - \S X \»(S, • ^ + ^o) - v lS ■ A S xx )dxdy 

= Dt JJ di^ylx, ^xlx) - ^v \lx\ 2 g - vig{Jlx, Vxlx)dxdy = (7.26) 

under suitable boundary conditions (where C denotes the coordinate domain of (x,y)). The 
vector field £(7) on S 2 , or equivalently the vector function £(S), in the non-polynomial 
constant of motion (7.23) has the geometrical meaning of a homothetic vector with respect 
to the metric-normalized volume form e g on S 2 , i.e. C^e g = e g . 

Theorem 7 is established as follows. In parts (i) and (ii), the derivation of the equations of 
motion (7.10) and Hamiltonians (7.12) for S(t,x) involves combining the evolution equation 
(6.3) for the frame vector E-i — S with the hierarchy (6.11) for the variable w by means of 
the identities 

u(E 2 - iE 3 ) = S x - iS A S x , (7.27) 
u y (E 2 - iE 3 ) = -(S A S(S y ) + iS(S y )), (7.28) 

in addition to 

Re(u/) =S x -f, lm(uf) — (S A S x ) ■ f, (7.29) 
Re(u y f)=S(Sy)-(SAf), lm(u y f)=S(Sy)-f (7.30) 
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holding for vectors / in R 3 , with components / = (E 2 + 1E3) • f, such that / • S = 0; here 
S is the recursion operator (4.24). Similarly, the derivation of the Hamiltonian structures 
(7.17) in part (iii) for n 7^ 1 relies on applying the previous identities to the bi-Hamiltonian 
structure (6.13) for the flow equations on u(t,x). The n = 1 case reduces to computing the 
Hamiltonian derivatives (7.21)-(7.22), which is carried out in appendices B and C. Finally, 
all of the corresponding results for j(t,x) are an immediate consequence of the geometric 
identifications (7.5)-(7.7) in addition to (3.41)-(3.42). 

VIII. GEOMETRIC FORMULATION 

There is a natural geometric formulation for the surface flows (6.1) corresponding to 
the 2+1 vector models (7.10) and 2+1 dynamical maps (7.11) in Theorem 7. We begin 
by writing down the the intrinsic and extrinsic surface geometry of f(x, y) in terms of the 
variables u,v,q,v±,qi appearing in the structure equations of the parallel framing for the 
non-stretching x-coordinate lines. 

Proposition 7. Let f(x,y) in IR 3 be a surface with a parallel framing (5.1)-(5.5) adapted 
to the x coordinate lines, satisfying the structure equations (5.16) and (5.17). Then, on the 
surface r(x,y), the infinitesimal arclength is given by the line element 

ds 2 = (dx + qidy) 2 + \q\ 2 dy 2 , (8.1) 

and the infinitesimal surface area is given by the area element 

dA=\q\dxAdy. (8.2) 

All other aspects of the intrinsic surface geometry can be derived from the line element 
(8.1). In particular, the 1st fundamental form (i.e. the surface metric tensor) is simply 
dxdx + 1q\dxdy + \q\ 2 dydy, from which the Gauss curvature can be directly computed in 
terms of the x,y coordinates [12]. 

The extrinsic surface geometry can be determined through the surface normal vector 

ft = f x A f y = lm(q(E 2 + iE 3 )), \n\ = \q\, (8.3) 

as given by the expression (5.13) for the surface tangent vectors f x and r y in terms of 
the parallel frame along the ^-coordinate lines. This normal vector (8.3) depends on a 
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choice of the transverse coordinate y due to its normalization factor \q\. To proceed, we use 
the following natural geometric framing [12] that is defined entirely by the non-stretching 
direction on the surface and the orthogonal direction of the surface unit-normal in M 3 : 



e 11 = f x = Ei, 

e L = Inl^n = lm(e^(E 2 + iE 3 )), 

*e" = Inl" 1 ^ Af x = Re(e # (£ 2 + iE 3 )), 



(8.4) 
(8.5) 
(8.6) 



where 



i> = arg(g). (8.7) 

(Here * denotes the Hodge dual acting in the tangent plane at each point on the surface.) 
Note this frame (8.4)-(8.6) differs from a parallel frame by a rotation through the angle 
(8.7) applied to the frame vectors in the normal plane relative to the x-coordinate lines, 
with e" and *e" being an orthogonal pair of unit-tangent vectors on the surface f(x,y), and 
e 1 - being a unit-normal vector for the surface. 

The Frenet equations of the frame e", *e", e 1 - directly encode the extrinsic geometry of the 
surface f(x, y). In matrix notation, the Frenet equations with respect to the x, y coordinates 
are given by 

/ n „, „\ / .1 \ / .1 \ / 



*e 
i 

e 



a (3 
—a 7 
-P - 7 



*e 



\ 6 Jy 



where 



a 



*e" 



Ref-ue 



n p 
—fi a 
-p -a } 

f3 = e ± -el = Im(Me^), 



*e 



7 = e ± ■ *el = \q\ 1 lm(q x e ii ') = -if> x 



(8.9) 
(8.10) 



and 



fj, = *e 11 • ej = Re(ve 1 ^) = q x a + \q\ x , p = e x • ej = lm(ve^) = qrf - \q\ip x , (8.11) 
a = e ± ■ *ej = \q\' l lra((q y + iv iq )e^) =v 1 -^ y (8.12) 

are obtained through the relations (5.8), (5.9), (5.14) together with the structure equations 
(5.16), (5.17). With respect to the x,y coordinates on the surface, the scalars a and p are 
known as the geodesic curvatures; (3 and p as the normal curvatures; 7 and a as the relative 
torsions [12]. 
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Proposition 8. For a surface r(x,y) with a parallel framing (5.1)-(5.5) adapted to the x 
coordinate lines, satisfying the structure equations (5.16) and (5.17), the 2nd fundamental 
form is given by 

II = (f x dx + fydy) ■ (e^dx + e y dy) = —((3dxdx + 2pdxdy + (a — qi^)\q\dydy). (8.13) 

The components of II with respect to the surface tangent frame e",*e" yield the extrinsic 
curvature scalars 

k u = e" • J D| ! e ± = -p, (8.14) 

ki2 = e" • D^e 1 - = k 2 i = *e" • = -7, (8.15) 

k 22 = *e 11 ■ D^e L = -|g| _1 (o- - ^7), (8.16) 

where D\\ = e^\D = D x and D*\\ = *e^\D = \q\~ 1 (D y — q\D x ) denote the projections of the 
total exterior derivative D on the surface in the directions tangential and orthogonal to the 
x-coordinate lines. 



All aspects of the extrinsic surface geometry of f(x, y) can be determined from the ex- 

k n ki2 
&21 ^22 



k\\ K\2 

trinsic curvature matrix I . In particular, the mean curvature of the surface [12] 



H = {k u + k 22 )/2 (8.17) 

is given by the normalized trace of the extrinsic curvature matrix. 

Now, in terms of the frame vectors (8.4)-(8.6), any surface flow (6.1)-(6.2) in which the 
x-coordinate lines are non-stretching can be written in the form 

f t = ae 11 + b*e 11 + ce ± , \f x \ = 1 (8.18) 

with 

a = D- 1 Re(uh), b = Re(e^h), c = Im(e^/i). (8.19) 

Through Proposition 5, we then obtain a hierarchy of flows on r corresponding to the bi- 
Hamiltonian 2+1 flows on u in Theorem 6, as given by 

6 W + ic (n) = _ e <^n-i( Ui/ ) > a (») = -D^Re^"" 1 (u y )) (8.20) 
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in terms of the recursion operator 1Z = i(D x + uD x 1 Re(uC)). Moreover, these coefficients 
(8.20) have a geometrical formulation derived from the operator identity 

e^TZe-^ = iD x +ip x + iue^ D^Reiue^ C) = iD x - 1 -(/3- ia)D- l Re((a + i/3)C) (8.21) 

combined with 

e^u y = (e ilp u) y - iip y e ilp u = a y + /% + i(P y - atp y ). (8.22) 
This leads to the following geometric counterpart of Theorem 7. 

Theorem 8. The integrable 2+1 vector models (7.10) and integrable 2+1 dynamical maps 
(7.11) correspond to a hierarchy of surface flows in M. 3 , 

f t = a^-^e 11 +6 ( "- 1) *e 11 +c (n " 1) e ± , \f x \ = 1, n = l,2,... (8.23) 
with geometric coefficients 

6 W + ic (n) = - V ^\ ay + ^ y + i(p y _ Qr0y)) , a (n) = D~\ab^ + /3c (ri) ), (8.24) 
given m terms of the recursion operator 

<p = iD x - 1 -((3 - ia)D- l Re((a + i(3)C). (8.25) 

Here a and (3 are the geodesic curvature and normal curvature of the non-stretching x- 
coordinate lines on the surface r(t,x,y), 7 is the relative torsion of these lines, and ip = 
-j-fdx. 

Remark: The bottom flow in this hierarchy can be written in the alternative form 

&(°) = - p , c (°) = p, a® = D~\pn - ap) = v x (8.26) 
by means of the relation 

V-\ay + + i(P v - a^jy)) = p-ip (8.27) 

expressed in terms of the geodesic curvature p and normal curvature p of the y-coordinate 
lines on the surface r(t,x,y). This relation (8.27) is obtained from —ie l ^v = V~ x {e % ^u^) 
which is a straightforward consequence of Lemma 1. 
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Geometric properties of the surface flows in the hierarchy (8.23) can be straightforwardly 
derived from the results in Propositions 7 and 8 combined with the explicit evolution equa- 
tions for the variables q and v (as determined by Lemma 1). In terms of the surface flow 
equation (8.18)-(8.19) and the operator identity (8.21), these evolution equations are given 
by 

i Vt = WlV + e -^(D y + icr)V(b + ic), (8.28) 
qt = aq x - iqw 1 + e~ ilp (D y - q x D x + %{a - qi^)){b + ic), (8.29) 

with 

w lx = Re((a-i/3)V(b + ic)), q lx = \q\a = -Re((a - 0)V(p - in)), (8.30) 

from which we obtain the evolution of the geodesic curvature, normal curvature, and relative 
torsion of the x-coordinate lines: 

a t = D x {aa) + Re{V 1 2 {b + ic)) + (31m{V 2 {b + ic)), (8.31) 

fit = D x (af3) + Im(XV(& + ic)) - alm(V 2 (b + ic)), (8.32) 

7i = D x (ari + Im(P 2 (6 + ic))) - Re((/3 + ia)V ± (b + ic)), (8.33) 

where 

V l = D\\ - ik 12 , V 2 = D m \\ - ik 22 (8.34) 

are a pair of geometric derivative operators associated with the x-coordinate lines on the 
surface (cf. Proposition 8). Similarly, we find the area element of the surface has the 
evolution 

(dA) t = a(dA) x + Re{V 2 {b + ic))dA = edA + C x dA (8.35) 

which consists of an infinitesimal change due to the tangential part of the surface flow X = 
ae" + 6*e'l plus a multiplicative expansion/contraction factor e = 2Hc related to the mean 
curvature (8.17) of the surface through the normal part of the surface flow C e ±dA = 2HdA. 
These developments now lead to the following geometric results. 

Theorem 9. Under each flow n = 1,2, .. . in the hierarchy (8.23), the surface r(t,x,y) is 
non-stretching the x-coordinate direction while stretching in all transverse directions, such 
that surface area locally expands /contracts by the dynamical factor 

e (n) = 2Hlm{V n ' 1 {p-ifi)) 
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where H is the mean curvature (8.17) of the surface. The geodesic curvature a, normal 
curvature (3, and relative torsion 7 = — ip x of the x-coordinate lines in each surface flow 
(8.23) satisfy the integrable system of evolution equations 

a t + i/3 t = D x {a^- 1) {a + ip))-V 1 2 V n - 1 {p-iii) - (f3 - ia)Im(I? 2 P n - 1 (p - ip)) (8.36) 
i> t = a^V* + a (n) + lm(V 2 V n - l (p - ip)) (8.37) 

in terms of the recursion operator (8.25) and the pair of geometric operators (8.34), w ^h 
the geodesic curvature p and normal curvature p of the y-coordinate lines given by (8.27), 
and where 

a W = - J D~ 1 Re((a - i/3)V n (p - ip)) (8.38) 

yields (modulo total x,y -derivatives) a set of non-trivial constants of motion JJ c a^dxdy = 
IJJ C a(3 y - f3a y - ip y (a 2 + (3 2 )dxdy, J J c a^dxdy = \$J C a x a y + (3 x (3 y + ifj x ip y (a 2 + (3 2 ) + 
i/j x (a y f3 — f3 y a) + ip y (a x f3 — (3 x a)dxdy, etc. for the system (8.36)-(8.37). (Here C denotes 
the coordinate domain of (x,y)). 

We conclude by pointing out that the surface flows (8.23) in Theorem 8 and the corre- 
sponding integrable systems (8.36)-(8.37) in Theorem 9 provide a geometric realization for 
the 2+1 vector models (7.10) and 2+1 dynamical maps (7.11) in Theorem 7. 

Ex. 1: The +1 surface flow is given by 

n = v x r x + pf x A n + ph, \f x \ = 1, (8.39) 

with V\ x = (5p — ap, where n = ^/\r y \ 2 — (f x ■ r y ) 2 f x A f y denotes the surface unit-normal. 
This flow is a geometric realization of the 2+1 Heisenberg model (7.1), corresponding to the 
new integrable 2+1 generalization of the Schrodinger map (7.8). 

Ex. 2: Similarly, a geometric realization of the 2+1 mKdV vector model (7.3) and the 
corresponding mKdV map (7.9) is provided by the +2 surface flow 

- f t = v f x + (a y + i> y P)r x A n + (/3 y - ip y a)n, \f x \ = 1, (8.40) 

with u 0x = aa y + (3f3 y . 

The corresponding geometric integrable systems on the geodesic curvature a, normal 
curvature j3, and relative torsion 7 = — ip x of the non-stretching x-coordinate lines in these 
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surface flows (8.39) and (8.40) are respectively given by 

a t + ip t = v x (a x + i/3 x ) - V x 2 (p - ip) + (a + i/3)(p/3 - pa + ilm(V 2 (p - ip))) (8.41) 
i> % = v^ x - u + lm(V 2 (p - ip)) (8.42) 

and 

a t + if3 t = -u (a x + if3 x ) - V 1 2 (a v + + i(P y - OLip y )) 

-(a + i/3) (i(a 2 + f3 2 ) y - ilm(V 2 (a y + /% + i(/3 y - a^ y )))) (8.43) 
ipt = -iso*Px + V2 + Im(D 2 (a y + f3^ y + i(f3y-a^ y ))) (8.44) 

with v 2x = af3 xy - f3a xy - ip y {aa x + f3f3 x ) - ip x {aa y + (3f3 y ) - ip xy {a 2 + f3 2 ). 

IX. CONCLUDING REMARKS 

There are some directions in which to extend the geometrical correspondence among 
integrable vector models, Hamiltonian curve and surface flows, and bi-Hamiltonian soliton 
equations presented in this paper. 

First, it would be of interest to generalize this correspondence to integrable models for 
spin vectors S in Euclidean spaces M. N for iV > 3. In particular, the Heisenberg model 
S t = S AS XX , S-S — 1, is well known to have a natural generalization where the vector wedge 
product and dot product in IR 3 are replaced by a Lie bracket [ , ] and (negative definite) 
Killing form ( , ) of any non-abelian semisimple Lie algebra on M N , i.e. St = [S, S xx ] , 
— (S,S) = 1. All of our results in this paper have a direct extension to such a Lie algebra 
setting by applying the methods of Ref.[5] to non-stretching curve flows in semisimple Lie 
algebras viewed as flat Klein geometries. This will lead to a large class of integrable 2+1 
generalizations of the Heisenberg model for spin vectors S in semisimple Lie algebras. 

Second, an interesting open problem is to find a similar geometric derivation for non- 
isotropic spin vector models in IR 3 as described by the Landau-Lifshitz equation St — S A 
(S X x + JS), S ■ S = 1, where J = diag(ji, j 2 , j 3 ) is a constant matrix which measures the 
deviation from isotropy. This equation has two compatible Hamiltonian structures [6], one 
of which uses the same Hamiltonian operator SA that arises in the isotropic case (i.e. in 
the Heisenberg model). The second Hamiltonian operator, however, involves the anisotropy 
matrix J, which cannot be derived from the frame structure equations for non-stretching 
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curve flows in Euclidean space. This suggests a non-Euclidean geometric setting will be 
needed instead. 
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Appendix A: Proof of Theorem 2 

Let f(x) be a space curve with x as the arclength, i.e. \f(x)\ = 1, and let T, N, B be its 
Frenet frame (2.10). To prove Theorem 2, we will enumerate the Euclidean invariants of 
f(x). 

Firstly, since the curvature k = T x -N = \T X \ and torsion r = N x -B = \T X \~ 2 T XX - (T AT X ) 
are invariantly defined in terms of the unit tangent vector T — f x along r(x), so are all of 
their x derivatives. This establishes part (iii) of the theorem. 

Secondly, these (differential) invariants generate all possible scalar expressions formed out 
of T, T x , T xx , ... by dot products and wedge products, as shown from a recursive application 
of the Frenet equations (2.14)-(2.15). Specifically, 

T x = kN, (Al) 
T xx = —k 2 T + k x N + ktB, (A2) 
T xxx = -3kk x T + (k xx - k 3 - kt 2 )N + (2k x t + t x k)B, (A3) 
etc. 

yields 

T x ■ T x = —T xx • T — k 2 , (A4) 
1 1 

T xx • T x —(T x ■ T x ) x ~^T XXX ■ T nn X j (A5) 

T xxx • T x —(T x ■ T x \ xx T xx ■ T xx ~^T XXXX ■ T ^^xx ^ k t , (A6) 
etc. 
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and 



T xx ■ (T A T x ) = k 2 t, 

T xxx • (T A T x ) = (T xx • (T A T x )) x = k 2 t x + 2tkk x , 

T xxx ■ (T A T xx ) = kt(k xx k kt ) + k x (2k x t + t x k), 



(A7) 
(A8) 



(A9) 



etc. 



which thus establishes parts (i) and (iv) of the theorem. 

Finally, on the other hand, since a parallel frame along f(x) is unique up to a rigid (x- 
independent) rotation on the normal vectors, the corresponding components of the principal 
normal T x of r{x) given by u = T x ■ E L = ne 10 are invariantly defined only up to rotations 
9 — >■ 9 + 0, with =const., on E 1 = (N + iB)e ld (and likewise for the components of T xx , 
T xxx , etc.). Such U (1) rotations comprise all transformations preserving the parallel property 
(2.20)-(2.21) of this framing. Consequently, the actual invariants of r(x) will correspond to 
[/(l)-invariants formed out of u, u, u x , u x , . . . via the relations k = \u\, r = (axgu) x . This 
establishes part (ii) of the theorem. 

Appendix B: Hamiltonian structure of the 1 + 1 and 2+1 Schrodinger maps 

We will first verify the second Hamiltonian for the Schrodinger map equation (4.54) and 
its 2+1 generalization (7.20). The following preliminaries concerning the tangent space 
structure of S 2 will be needed. Here u,v,w will be any triple of tangent vectors. 

(1) The metric tensor g, complex structure tensor J, and metric-normalized volume form 
t g on S 2 satisfy the identities 



(2) In local coordinates on S 2 , the metric-compatible covariant derivative V (i.e. Riemannian 
connection) and the associated covariant divergence operator div g are given by 



g(u, Jv) = € g (u,v), 
€ g (u, v)w + € g (v, w)u + € g (w, U)V = 0. 



(Bl) 
(B2) 



V v u = d v u + F v u, 



(B3) 



div g u = divu + tr(r„), 



(B4) 
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where T denotes the Christoffel symbol [12] determined from g by the properties 

T v u = T u v, (d w g) [u, v) = g(v, T w u) + g(u, T w v) . (B5) 

(3) For an arbitrary variation £7 of the map 7 into S 2 , geometrically represented by a tangent 
vector field, the variation of g\ 7 and e 9 | 7 (as induced by their evaluation at 7) is given by 

6gh(u, v) = d Sj g(u, v) = g(u, T Sy v) + g(y, T Sj u), (B6) 
5e g \j(u,v) = (y/det g^ds^y/det g)e g (u, v) = tr(r 57 )e 9 (w, v). (B7) 

Now, consider the 1+1 Hamiltonian density (4.52) defined in terms of a vector field £(7) 
with covariantly-constant divergence, div s £(7) = I. Through the identity (Bl), this density 
can be written more conveniently as 

H^ = eMl),l*)- (B8) 

Its variation is given by 

= tr(r 57 )e g (e( 7 ),7,) + ^£(7), 7.) + ^(7),^)- (B9) 
Integration by parts on the third term in (B9) yields 

e 9 (£( 7 ), DJ&i) = D x {eMl), ' ^£(7), *y) - tr(r 7 Je fl (£( 7 ), *y)- (BIO) 
By combining the middle terms in (B9) and (BIO) via the identity (B2), we get 

tgidsytil), Ix) ~ e 9 (djj(7), h) = e 5 (^7, 7x)div^( 7 ). (Bll) 
Likewise, combining the first term in (B9) with the third term in (BIO), we obtain 

e fl (£( 7 ),7*)tr(r 57 ) - e g (£( 7 ),<*7)tr(r 7 J = e fl (5 7 , 7 *)tr(r €(7) ). (B12) 
Hence, modulo total x- derivatives, (Bll) and (B12) combine to give 

^(-1) = e fl (5 7 , 7x )div ff f (7) = g(6>y, J lx ) 

which yields the Hamiltonian derivative (4.53). 

Next, consider the 2+1 Hamiltonian densities (7.13) and (7.14). The previous derivation 
applies verbatim to the density (7.13), yielding its derivative (7.21). For the second density 
(7.14), we look at its two terms if — H 1 + H 2 separately: 

Hi = ^g(lx,ly), H 2 = 7^ie fl (£(7),7x), (B13) 
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with 

vix = e 9 (7x,7j/)- (B14) 
First, the variation of Hi is given by 

SH i = ^9{lx, D y 5~i) + ^g(DJ^, 7l/ ) + ^g(^ x , T s ^ y ) + ^(t v , r«5 7 7x) 

= ^(7x, v y< J 7 ) + lady, V X< J 7 ) (B15) 

through equations (B3) and (B5). Integration by parts on these terms yields 

girtx, V^ 7 ) = Dy{g{-f x , <S 7 )) - ^(fry, V^), ff( 7w , V x <& 7 ) = D x {g{-f y , <5 7 )) - #(5 7 , V x7y ) 

(B16) 

where, on scalar expressions, a covariant derivative reduces to an ordinary total derivative. 
Hence, modulo total x, y-derivatives, substitution of (B16) into (B15) yields 

5Hi = --g(6i, V xly + V ylx ) = g{5 1 , -V ylx ) (B17) 

after we use the commutativity identity V x ^ y = V y ^ x which is a consequence of the first 
property in (B5). 

Second, the variation of H 2 is given by 

5H 2 = leMl)^)8vi + l -ViD x {eMl)^l)) + ^(^7, 7^div 9 £( 7 ) (B18) 

where the last two terms come from (B10)-(B12). To evaluate the first term in (B18), we 
start with 

Sv ix = e g (DJ>y, 7y ) + e s ( 7x , D y 5>y) + tr(r 57 )e ff (7 x , 7y ) (B19) 
and use integration by parts to expand the first two terms, giving 

e g (DJ^, 7y ) = D x {e g {5-f, 7y )) - e 9 (5 7 , D x ^ y ) - tr(r 7 Je 3 (5 7 , 7y ), (B20) 
tgilx, D y 5-f) = D y {e g {-f x , 57)) - e g {D y -f x , <5 7 ) - tr(r 7y )e 9 ( 7:c , 57). (B21) 

Using the identity (B2), we note the Christoffel terms in (B19), (B20), (B21) combine to 
give 0, while the middle terms in (B20) and (B21) cancel due to 

D xly - D ylx = V xly - V ylx = 0. (B22) 

Hence, (B19) simplifies to a sum of total x, ^/-derivatives 

5v lx = D x (e g (5-f, 7w )) + D y (e g (-f x , 5 7 )). 
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As a result, the first term in the variation (B18) becomes 

\eMl)^)^i = ^7,7^(7), 7*) + ^/(Dye^d^eMl)^)- (B23) 
Integration by parts on the second term in (B23) yields 

^Dy^^^eMl)^) = \^MDZ X D v {eMi)^))- (B24) 
Now we use the relation 

DyCgitd), lx) ~ D x e g (C(-f)n y ) = e g (d ly £(l),lx) ~ ^g(d lx C(l)n y ) + ^(£(7), D ylx - D xly ) 

+e fl (e( 7 ),7.)tr(r 7 J - e fl (e(7),7„)tr(r 7 J 
= e s (7 y ,7 :c )(div^(7) +tr(r 5(7) )) = -uixdiv fl f (7) 

obtained via the identities (B2) and (B22). Thus, (B24) reduces to 

g^A/ta (7x, *7))e 9 (^(7), 7x) = )fg(lx, $l)tg(t(l),ly) + \ v itg(h, lx) (B25) 
which combines with the first term in (B23) by use of (B2) to give 

^(£(7), 7x)<^i = ^9(^(7), h)£g(lx, ly) + ^ie#7, 7x)- (B26) 
Substituting (B26) into the variation (B18), and using (B14), we get 

5H 2 = 2^es(C(7)^7) + 2 VlD *( e 9(Z(7), + Vi6 g (5-f,-f x ) = giSj^J^). (B27) 
Finally, we combine the separate variations (B27) and (B17) to obtain 

SH^ = SH, + 5H 2 = g{5 1 , -W ylx + ^ J lx ) (B28) 
which yields the Hamiltonian derivative (7.22). 



Appendix C: Hamiltonian structure of the 1+1 and 2+1 Heisenberg models 

We will next verify the second Hamiltonian for the 1+1 Heisenberg model (4.55), given 
by the density 

H ( ^=aS)-(SAS x ). (CI) 
Here £(S) is a vector function, defined in terms of the spin vector S, such that 

S-t(S)=Q (C2) 
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and 

dj-as) = i, (C3) 

where dg = ds — S(S ■ ds) is a component-wise gradient operator satisfying the properties 
(4.57). We note that, due to these properties, d$ has a well-defined action on any function 
of S with S ■ S — 1. To proceed, the following algebraic preliminaries will be needed. 

(1) A variation of S consists of an arbitrary vector SS _L S, i.e. S • SS — 0. 

(2) The variation of £(S) induced by SS is given by 

5t(S) = (5S.dj;)t(S). (C4) 
Similarly, the total x-derivative of ^(S*) is given by 

D^(S) = (S x -dknS)- (C5) 

(3) Since the subspace of vectors orthogonal to 5 in M 3 is two-dimensional, A S x lies 
in the one-dimensional perp space, so thus 

5S A S x = -{5S ■ {S A S X ))S. (C6) 

Now, the variation of the density (CI) is given by 

6H(-V = 6£(S) ■ (S A S x ) + ^(S) ■ (SS A S x ) + £(S) ■ (S A DJS). (C7) 

Integration by parts on the third term in (C7) yields 

aS) ■ (S A DJS) =D X (S- (SS A aS))) - aS) ■ (S x A SS) - SS ■ (D x aS) A S) (C8) 

with the middle terms in (C8) and (C7) each vanishing due to (C6). Hence, modulo total 
x- derivatives, (C7) reduces to 

6H(-» = SaS)-(SAS x )-SS-(D x aS)AS) = S-(S x A(SS-di)aS)-SS A(S x -dj)aS)) (C9) 

via (C4) and (C5). To simplify the terms in (C9), we first rewrite 

S X (SS ■ dj) - SS(S X ■ dj) = (SS A S x ) Adj = -SS ■ (S A S X )(S A dj) (CIO) 

by means of standard vector cross-product identities in addition to identity (C6). Thus, 
(C9) becomes 

SH ( - 1] = (SS ■ (S A S X ))S ■ (-(S A dj) A as)) (Cll) 
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and we again apply vector cross-product identities to rewrite the term 

(s a oj) a as) = dj(s ■ as)) - asm ■ s) - s(dj . as)). (cu) 

Then, since S ■ d$ = = S ■ we have 

S ■ (-(S A dj) A aS)) = dj ■ aS) = 1 (C13) 
whence (Cll) simplifies to 

SH^-V = SS ■ {S A S x ) 

yielding the Hamiltonian derivative (4.58). 

The above derivation carries over verbatim to also verify the first Hamiltonian (7.13) 
for the 2+1 Heisenberg model (7.19). To verify the second Hamiltonian (7.14), we will 
separately consider the two terms in the density H^ 1 ^ = Hi + H 2 given by 

H l= l -S x -S y , H 2 = ^ Vl aS)-(SAS x ) (C14) 

with 

v lx — S x • (S A S y ) — S • (S y A S x ). (C15) 
The following identity will be useful: 

S y A S x = v lx S (C16) 

holding similarly to (C6). 

First, the variation of Hi is simply 

5H X = ^S x ■ D y 5S + ±S V ■ DJS = -SS ■ S xy (C17) 

modulo total x, y- derivatives . Next, the variation of H 2 consists of the terms 

5H 2 = ^(S) ■ (S A S x )5 Vl + ^ Vl D x (6S ■ (^(S) A S)) + l - Vl 5S ■ (S A S x ) (C18) 

as obtained from (C8), (CIO), (C12), (C13). To evaluate the first term in (C18), we note 
the variation of (C15) is given by 

5 Vlx = 5S ■ (S y A S x ) + S ■ (S y A D X 5S) + S ■ (D y 5S A S x ) 
= D y (5S ■ (S x A S)) + D X (5S ■ (S A S y )) + 2v lx S ■ 5S 
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through (C16). Since the last term vanishes due to the orthogonality SS _L S, this yields 

5 Vl = 5S ■ (S A S y ) - D- 1 D y {SS ■ (S A S x )). 
Hence, the first term in (C18) becomes 

^(S)<SAS x )S Vl = ^5S-(SAS y )aS)-(SAS x )-^5S-(SAS x )D- 1 D y (aS)-(SAS x )) (C19) 
after integration by parts. We simplify the second term in (C19) by using the relations 
D y (t(S) ■ (S A S x )) - DMS) ■ (S A Sy)) 

= S ■ (S X A (Sy ■ d£)£(S) — Sy A (S X ■ dj)^S)) + 2£(S) ■ (Sy A S x ) 

where, similarly to (CIO) and (C12), 

S X (Sy ■ dj) ~ Sy(S X -dj) = S- (Sy A S X )(S A dj) 

yields 

S-(S X A (S y ■ dj)Z(S) — Sy A (S x ■ dj)Z(S)) = -v lx dj ■ ^(S), 

while 

^(S) ■ (Sy A S X ) — o 

holds due to (C16) and (C2). Thus, we have 

D?D y (t(S) ■ (S A S x )) = aS) ■ (S A S y ) - v l7 
whence (C19) simplifies to 

^(S) ■ (s a s x )5 Vl = Us ■ (s a s y ) as) ■ (s a s x ) - Us ■ (s a s x ) as) ■ (S A S y ) 

+U 1 5S-(SAS X ). (C20) 
By applying vector cross-product identities to the first two terms in (C20), we get 



^(5S A S) ■ Sy(aS) AS)-S X -^(5SAS)- SMS) A S) ■ S y 



^(5S A S) ■ ((S x A S y ) A (as) A S)) = ~(6S A S) ■ aS)v la 



via (C16). Hence, (C20) reduces to 



^(S) ■ (S A S x )5 Vl ee ±6S ■ (as) A S)v lx + ± Vl 6S ■ (S A S x ). (C21) 
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Finally, combining (C21) with the middle term in (C18), we get a total x-derivative, so thus 
(C18) becomes 

5H 2 = SS- ( Vl SAS x ). (C22) 

The separate variations (C22) and (C17) then combine to yield the Hamiltonian derivative 
(7.22). 
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